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The dielectric function and plasmon modes of a two-dimensional electron gas (2DEG) are studied 
in single- and double-quantum-well structures with Rashba spin-orbit interaction (RSOI) in the 
framework of the random-phase approximation. The RSOI splits each parabolic energy subband 
of a 2DEG into two nonparabolic spin branches and affects the electronic many-body correlation 
f*^ ■ and dielectric properties of the 2DEG. The influence of the RSOI on the 2DEG plasmon spectrum 

in single quantum wells appear mainly in three ways: 1) an overall frequency lowering due to the 
energy band deformation; 2) a weak frequency oscillation stemming from the spin-split energy band; 
and 3) an enhancement of the Landau damping as a result of the emerging of the inter-branch single- 
particle-excitation spectrum. In double quantum wells, the above effects are enhanced for the optic 
plasmon mode but diminished for the acoustic one. 
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INTRODUCTION 



'■^J , The growing interest in spintronics is leading to an extensive exploration of electronic properties of a two-dimensional 
£h ' electron gas (2DEG) in semiconductor heterostructures with spin-orbit interaction (SOI) thanks to its promising 
application in manipulating spins. p], Q As a relativistic effect of dynamic electrons moving in an electric field, 
1 the intrinsic SOI exists in bulk semiconductors with structural inversion asymmetry, Q while controllable SOI can 

be introduced by asymmetrically confining a 2DEG in semiconductor heterostructures to create an average electric 
t— I , field across the electron system. In general, the former results in the cubic Dresselhaus term of the SOI and the 
latter introduces the Rashba term for isotropic 2DEG and the anisotropic linear Dresselhaus term. 0, 0, |(| In some 
narrow-gap semiconductors, such as InGaAs, the Rashba spin-orbit interaction (RSOI) dominates and can be well 
controlled in the laboratory. Q In the past several years many authors addressed the RSOI effects on one-body 
ballistic transport but the study of the RSOI effects on many-body properties of 2DEG is relatively limited. As is 
well known, dielectric properties and plasmon spectrum can reveal fundamental many-body correlations in a 2DEG. 
In the presence of the RSOI they have been discussed in the literature. H , 0, IToL [ill lisj However, some authors 
concluded that the charge plasmon spectrum is not affected by the RSOI |T(J and others declared the splitting of the 
■ plasmon spectrum in the presence of RSOI. 0, 

From an experimental point of view, the double-quantum-well (DQW) structure is a good candidate for studying 
the electron correlation and collective modes of a 2DEG through transport measurements such as Coulomb drag. In 
a DQW system composed of two spatially separated quantum wells without tunnelling, the two 2DEG's, which are 
coupled with each other via Coulomb interaction, may oscillate collectively out of phase or in pha se. The former 
case corresponds to the lower frequency acoustic plasmon and the latter to the optic plasmon. |l3l IbH ITU ITri Il7j 
Both plasmon modes play an important role in the many-body properties and enhance significantly the Coulomb 
drag in bilayer systems. [Lsl IT^. I20L I2H ] In this paper we focus on the effects of the RSOI on the plasmon modes in 



single-quantum-well (SQW) and DQW systems. 

II. DIELECTRIC FUNCTION 

We consider two Coulomb-coupled 2DEG's without tunnelling in a DQW structure of a narrow gap semiconductor, 
InGaAs/InAlAs, where the RSOI due to confinement along z direction dominates. An effective mass m* = 0.05too 5 
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k x (10 5 cm 1 ) 

FIG. 1: Energy spectrum of a 2DEG with RSOI. The two dotted ellipses at the Fermi energy Ef illustrate the rotational 
symmetry in x-y plane of the two Fermi circles of the + and — spin branches. The arrows tangent to the Fermi circles indicate 
the spin directions. 

with mo the free electron mass, and the background dielectric constant = 10.8 are assumed. In each well the 
one-electron Hamiltonian of the 2DEG, in the x — y plane, is expressed as 

The parameter a measures the RSOI strength and is proportional to the average interface electric field, V(z) is the 
electric confinement potential along the z direction, a — {cr x ,(jy,cr z ) denotes the spin Pauli matrices, and p is the 
momentum operator. The eigenstates of the Schrodinger equation H^(r) = E^(r) are readily obtained as 




with k = (k x , k y ), r = (x, y), and e 4 ^ k = cr(k y — ik x )/k; k = yk 2 + k 2 is the radial wavevector and a = ±1 denotes 

the electron spin in the + (upper) or — (lower) branch. Tp{z) is the wavefunction along the z direction. 0£ is the angle 
between the x-axis and the direction of the spin polarization, which is perpendicular to the electronic wavevector in 
the 2D plane as shown by the arrows along the Fermi circles in FigQ] 
The eigenvalues corresponding to the eigenstates (k, a) are 

E a M = ^ + aak. (3) 

This rotationally symmetric electron energy spectrum, E° vs k, is illustrated in Fig^for a strength a = 5ao and an 
electron density n e — 1.5no, with ao = 10~ n eVm and no — 10 n cm~ 2 . The Fermi wavevector of the spin branch a 
is kp = ^J2im e — k%, — ak a with k a = am* /H 2 . 

Eqs.(l)-(3) apply to each well when the inter-well coupling is neglected and a well index, i or j, will be introduced in 
the following when the coupling is taken into account. With the help of Eq.©, we obtain directly the bare Coulomb 
interaction between an electron (k, a) in well i and another one (p,o"i) in well j. As shown diagrammatically in 
FigEt a )> its Fourier transformation reads: 

<J, £,p 1,<T5 (q) = (k, k + qK (q)g^ (p, p - q), (4) 

where Vij(q) = J dre lqr / dz J dz'\tpi(z)\ 2 \ipj(z')\ 2 e 2 / (ineoei^/r 2 + (z — z 1 ) 2 ) is the bare Coulomb interaction between 
electrons without spin.eo is the vacuum dielectric constant and tjji(z) is the wavefunction along the z direction in well 
i, i = 1 or 2. As illustrated in Fig|5fb) , the spin introduces a nontrivial vertex g a ^ (k,k + q) = [1 + <ja'e l ^ k ^' t ' k -^]/2 
if the final spin polarization of electron differs from the initial one during a Coulomb collision. 
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FIG. 2: Diagrammatic representation of (a) the bare Coulomb potential v i (q), (b) the spin vertex g^ ,<T (k, k + q), and 

(c) the RPA Dyson's equation for the screened Coulomb potential of a 2DEG t^yj, ^ ri ' x (q, w). 



The technique for multicomponent system 0, is used to derive the dynamically screened Coulomb potential in 
the random phase approximation (RPA). The following self-consistent equation is found, as illustrated diagrammati- 
cally FigHc), 



<T2,cr 2 
m.s 



with 

U ^ (q ' U) = /*, + J5-(k)-i5"'(k + q)+W (6) 
and /i the Fermi distribution function in well i. 

! ! I I ! 

Note ET^-'h p 1,<T 1 (Qj w ) — sr* 7 k + G L)U i j(q,uj)g° 1 '' Tl (p, p — q), so the above equation reduces to 

Uij(q,w) =Vij(q) + ^2 Wi m (q).9,T^ CT2 (s,s-q)^^ 2 ; s ' T2 (q,w)5m 2 '' T2 (s-q,s)L7 m: ,■(q,w). (7) 

(J2,<?' 2 

m,s 

The elements of the corresponding 2x2 RPA dielectric matrix e, defined by v%j(q) — J2 m e im(q,u)U m j(q,Lj), reads 

eij(q,u) = 5ij - v i:j (q) ^ ^ (1 + aa 1 k ^ )nj;^ (q.u;) (8) 



CT,<T 

:/.k 



with being the angle between k and q. 



III. PLASMON DISPERSION 



The inverse of the dielectric matrix, e -1 describes the many-body response to the two-body bare Coulomb interaction 
Vij. At the poles of e _1 , i.e. for a set value of (q, w) a particular element of e _1 becomes singular, a density fluctuation 
of the 2DEG can lead to a density oscillation in the system or the onset of a collective mode. It is well known that the 
zeros of the real part of the determinant of e describes the plasmon dispersion and the imaginary part the damping 
of plasmon. In this paper, we discuss the plasmon spectrum by solving the equation Re\e(q,u>)\ = for the 2DEG in 
a SQW and a DQW in the presence of RSOI. 
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A. Plasmons in a single well 



If the coupling between 2DES's in different quantum wells is negligible, each well is independent and can be treated 
as a SQW. The well indices i and j drop and the dielectric matrix reduces to a scalar function 

e(q,u) = l-v(q)U(q 1 w) (9) 

with the nointeracting correlation function II(q, 10) = J2a a' n <T,cr (q, to) — ^ 5 ff , k [l + aa'(k + qcos9)/\k + 

q|]n^' CT (q, oS)/2, Although the RSOI splits the energy bands of different spins, the many-body Coulomb interaction 
between spin branches is described by a scalar dielectric function in the RPA rather than a matrix as in confinement 
induced multisubband systems. [22| The reason is that the RSOI does not change the wavefunction along z direction 
and the spin vertex is separable from the no-spin Coulomb interaction as shown in Eq. Q . Eq. coincides with the 
dielectric function obtained in Ref. and fill but not with Eq. (2) in Ref. IT2I 

In the absence of SOI, the zero-temperature plasmon dispersion relation, frequency u>o vs wavevector q, is approxi- 
mately expressed as 

^ = 4(nf/(A*- 1) ^ + q2{A " 1)] (10) 



with A = 1 + irh /[m*v(q)} and the Fermi wavevector fc^ = y/2i:n e . In the long wavelength limit q <C /cf, the 
above equation reduces to the well-known zero-temperature 2D plasmon frequency uj p — [n e e 2 q/(2eoCi , m*)) 1 / 2 if the 
thickness of the quantum well is negligible. 

In the presence of RSOI and q <SC k F , the intra-branch correlation function is given approximately by 

IT< CT (q, u) = (1 + ak a /k F )q 2 {k F ) 2 /{2nm*uj 2 ). (11) 

The second term is a result of the deviation of the energy band from parabola. Since k F > k F , this term leads to an 
extra negative term for the total intra-branch correlation II ++ + II and lowers the frequency of the plasmon. The 
inter-branch correlation is negligible in the long wavelength-limit because electron spins of the same wavevector in 
different branches are opposite. The long-wavelength 2D plasmon frequency with RSOI then also has a very simple 
dependence on the system parameters, 

W ° = Wj ,[l-A£/(47rne)]. (12) 

In the general case the plasmon dispersion relation is obtained numerically by searching the zeros of the real part 
of Eq. . For the sake of simplicity, only electrons in the lowest confined subband are considered and the electron 
wavefunction along the z direction is approximated by that of an infinite highly square-well potential with width b = 
23A, which is valid if the quantum well is not too wide. The bare Coulomb potential is then v(q) = e 2 F(qb)/(2eoeiq) 

withF(0 = {3t + 87r 2 /t~327r\l~ei)/ie(e+K 2 )}}/(e+K 2 )- 

In Fig. Ola) we plot the plasmon spectrum and the single continuum of a 2DEG with RSOI in an InGaAs/InAlAs 
SQW of typical parameters appeared in recent experiments. The plasmon spectrum is very close to that without 
SOI and it becomes difficult to distinguish one from the other. On the other hand, compared to the single continuum 
of the 2DEG without SOI, the single continuum with RSOI has an additional inter-branch band through which the 
plasmon spectrum passes. The real part of the dielectric function at q = 0.2 x 10~ 5 cm -1 is plotted in Fig. 01b), 
whose zero near u) = 7 meV gives the plasmon energy and zero near u) = 1 meV the edge of the intra-branch single 
continuum. Because the spin of an electron in the upper branch is exactly opposite to that in the lower branch 
at the same wavevector q, the vertical inter-branch transition is not allowed under the spin-independent Coulomb 
interaction, there is no zero corresponding to the edge of the inter-branch single continuum. Non- vertical inter-branch 
transitions, where the final electron wavevectors are different from the initial ones, are possible but very weak and 
have negligible effect to the plasmon spectrum for &p 3> k a or high electron densities as indicated in Eq. (|12fl and 
illustrated in Fig. [3] Nevertheless, since the plasmon spectrum passes through the inter-branch single continuum as 
shown in Fig. |3f a), the weak inter- branch transitions increase the Landau damping of plasmons and lead to a nonzero 
imaginary part of the dielectric function in the corresponding energy range. In Fig. OJc) the nonzero value of Im[e] 
in the energy range (6, 8.3) meV is attributed to the inter-branch transitions and is much smaller than the value of 
Im[e] inside the intra-branch single continuum in the energy range (0, 1) meV. 

The effect of RSOI on the 2D plasmon spectrum illustrated above coincides with that presented in Ref. 0, where 
the spectrum is obtained by the density-functional formalism and found to be independent of the RSOI. However, 
the result disagrees with the one obtained in Ref. |9j, where the optical conductivity is used to evaluate the plasmon 
spectrum and one extra mode is observed at the edges of the inter-branch single continuum. The disagreement 




FIG. 3: (a) Plasmon spectrum of a 2DEG of density n e = 20no in a SQW with SOI of strength a = Qo (solid curve) and in 
a SQW without SOI (filled squares). Inter- (darker shaded) and intra- (light shaded) branch electron-hole continuums in the 
presence of SOI are also shown. The vertical dashed line indicates the cross-section of the spectra at q — 0.2 x 10 -5 cm -1 , (b) 
The real part and (c) the imaginary part of the dielectric function e vs energy to at q — 0.2 x 10 -5 cm -1 . Limited by the space, 
the segments of the curves for Re[e] < -40 in the energy range (1,3) meV, Re[e] > 40 in (0,1), <Im[e] < 0.01 in (2.2,6.2) 
and (8, 10], and 10 <Im[e] in (0, 1.1) are not shown. The temperature is T — 2 K. 




qClOW) q(10 5 cm _1 ) 

FIG. 4: Plasmon spectrum of a 2DEG without SOI (dotted curve) and with RSOI (solid curve) in a SQW (a) and in a DQW 
(b) . The light (darker) shaded indicates the inter-branch (intra-branch) electron-hole continuum in the presence of RSOI . 

originates from the different approaches employed. In this paper, the dynamically screened Coulomb interaction is 
evaluated to find the plasmon modes and the vertical transitions between the two SOI branches are not allowed under 
the bare Coulomb interaction. As a result, the plasmon spectrum is not greatly affected by the SOI splitting. Similarly, 
the independence of the plasmon spectrum of the RSOI found in Ref. [I(J is a result of the fact that the charge density 
matrix elements vanish for vertical transitions. However, the nonvanishing matrix elements of the velocity operator 
v = VpH corresponding to vertical transitions introduces a new singularity to the optical conductivity. Consequently, 
a new plasmon mode is predicted in Ref. by an analytic treatment valid in the collision-free limit and involving 
order-of-magnitude estimates, see Eqs. (13) and (14) in Ref. Since the validity of this formalism in the presence of 
RSOI is ambiguous by itself, a direct judgement on it is beyond the range of this paper. Nevertheless, experimental 
observations of the plasmon spectrum of electrons in a quantum well with typical parameters as used in Fig. [21 is 
possible to clarify the controversy. 

In FigQJa), we plot the plasmon dispersion relation (solid curve) as well as the single particle excitation spectrum 
due to intra-branch (dark shadow area) and inter-branch (light grey area) transition for a 2DEG of electron density 
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FIG. 5: (a) Frequency ratio u)/uiq of a SQW plasmon vs the RSOI strength a for different wavevectors q = (solid curve), 10 5 
cm -1 (dotted curve), 2 x 10° cm -1 (dashed curve), and 3 x 10 5 cm^ 1 (dot-dashed curve). In the inset II _ ' + vs a for q = 0.1 
and u> = 3meV is shown, (b) ui/u>o vs q with values of a in unit of ao shown beside the curves. The filled (empty) circles show 
the resonant positions lu = 2akp (2akp). The electron density is n e = 1.5no in both panels. 

n e = 1.5no and RSOI strength a = 5ao- For the sake of comparison, the plasmon dispersion relation in the absence of 
RSOI is also shown (dotted curve). We observe an overall lowering of the plasmon frequency in the presence of RSOI, 
which results from the RSOI modification of the band structure. Different from the case without SOI where no overlap 
between the plasmon and the single-excitation spectra happens in the energy range shown in Fig^Ja), the plasmon 
spectrum with RSOI has a wider range of overlap with the interbranch single-excitation spectrum introduced by the 
RSOI spin split. Though numerical result shows that the imaginary part of the dielectric function is thousand times 
smaller in the inter-branch single-excitation regime than that in the range of intra-branch regime, this energy overlap 
can enhance greatly the Landau damping of 2D plasmon with RSOI. Near q — 2 x 10 5 cm -1 where uu{q) ~ 2akp, we 
notice also that a shoulder appears in the solid dispersion curve in Fig0Ja). Its explanation is given in the following. 

To have a closer look at the plasmon spectrum shift of RSOI, we define the frequency ratio lu/uuq as the plasmon 
frequency in the presence of RSOI divided by that in the absence of RSOI when other parameters are kept the same. 
In FiglSJa), the frequency ratio lu/luq is plotted as a function of the RSOI strength a for various wavevectors. In the 
long- wavelength limit q = (solid curve), only intra-branch transitions occour and lu/luq vs a follows Ea. (|12|) . For 
finite q (dotted and dash-dotted curves), inter- branch transitions also affect the plasmon spectrum especially near the 
inter- branch resonant energies at the Fermi wavevector of the — branch, Lu(q) ~ 2akp, where an increase of plasmon 
frequency may appear, and of the + branch, Lu(q) ~ 2akp, where a plasmon frequency decreases and a spectrum 
shoulder appears as shown in Fig. IHa). This behavior can be understood by checking the inter-branch correlation 
term in Eq. (JSJ. For positive lu, II + ' _ (q, lu) is negligible and II~' + (q, lu) is approximately expressed as 

ft-' + (q, u) = l<T' + (k, k - q)| 2 n k <+(q, w) cx q 2 k a /[k F (tkU ~ 2ak F )], (13) 
k 

with kp < kp < kp. A numerical result of II k ' + (q, u>) as a function of a is shown in the inset of Fig[5{a), which 
present smoothed peaks compared with that given by Ea. i|13|) . For strong RSOI, i.e. 2akp > u>, H~' + is positive and 
leads to an increase of plasmon frequency. If this positive inter-branch correlation is larger than the RSOI negative 
extra term of the intra-branch correlation in Ea. l|ll[l . there will be a RSOI enhancement of the plasmon frequency. 
This enhancement is observed in FigOJa) near a — 2ag with u/u>o > 1. For 2akp < u), Tl~ ,+ becomes negative and 
results in a further lowering of plasmon frequency. At fixed q, the frequency lowering of the inter-branch correlation 
appears generally stronger than the frequency enhancement effect because it happens at larger values of a and II~' + 
is proportional to a. In Fig. f5Ta). the curve of q = 10 5 cm -1 has two dips near a = 3«o and 9«o respectively because 
there are resonant energies lu ~ — E~ + = 2akp here. 

Kjp fx, 

Since the plasmon frequency ui in Ea. p3|) depends on the wavevector q, the frequency ratio lo/loq vs q shows similar 
frequency enhancement and lowering as illustrated in Fig|S{b). The effect of RSOI via intra-branch correlation is 
reflected in the value of lo/loq at q — and the contribution via inter-branch correlation is observed in the dependence 
on q. There is a lowering of the frequency ratio at low q and a frequency enhancement for high q. Each uj/loq curve 
reaches its minimum (maximum) at q about co(q) = 2akp (2akp). Among the curves in Fig. EJb), the effect of this 
inter-branch correlation is maximum for the curve a = 7ao whose minimum happens at the maximum value of q. 
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FIG. 6: Frequency ratios u>/u)q in a DQW as functions of a and q for the optic plasmon in (a) and (b), and the acoustic plasmon 
in (c) and (d). In panels (a) and (c), uj/ujo vs q for q = (solid), 10 J cm -1 (dotted), 2 x 10 5 cm -1 (dashed), and 3 x 10 5 cm" 1 
(dot-dashed) are shown. In panels (b) and (d) oj/wo vs q for (from up to down) a/ao = 1,2.5,3.5,5,7,8,9 are shown. The 
filled (empty) circles in (b) show the resonant positions as in Figl^Jb). The same electron density as in Figgis used. 

B. Plasmons in a double well 

In the DQW case, we present the results of a structure composed of two identical infinitely high quantum wells of 
width b = 23A separated by a distance d = 175A. The intra-well Coulomb interaction is vn(q) — v 2 2(q) = v(q) and 
the inter-well one v 12 {q) = v 2 i{q) = v d (q) = e 2 e - c " 1 G{qb)/(2eoe l q) with G(f) = 16tt 4 (1 - e«) e -V[£ 2 (£ 2 + ^?\- The 
plasmon spectrum is given by the zeros of the determinant of the dielectric matrix. From Eq.JSJ), the determinant is 
expressed as 

\e\ = {i-vii 1 ){i-vti. 2 )-v 2 d n 1 n 2 (14) 

with H the noninteracting correlation function in well i. In general, two plasmon modes, the optic mode of higher 
energy and the acoustic mode, are obtained out of the single particle excitation spectrum. If the electron densities in 
both wells are the same, Eq. (|10(l is still valid in the absence of RSOI for the zero-temperature plasmon dispersion by 
replace A with A op = v + Vd for the optic mode and A ac = v — Vd for the acoustic mode. In the long- wavelength limit, 
the frequency of the optic mode is about u>° p = \[2ui v and that of the acoustic mode ujp C = (n e e 2 d/2eoeim*) 1 / 2 q . 

In the numerical calculation, the plasmon dispersion is obtained by obtaining the zeros of the real part of the 
determinant of the dielectric matrix Re\e(q,uj)\. A typical plasmon dispersion of 2DEG in a DQW with RSOI is 
illustrated in FigQJb). As in the SQW, the dispersion curve of each plasmon mode is generally lowered by RSOI 
via intra-branch correlations and has a shoulder at the resonant energy 2akp of inter-branch transitions for electrons 
near the inner Fermi circle. Compared with the shoulder of the SQW plasmon spectrum in Fig0Ja), the shoulder 
for the optic mode happens at a lower wavevector by a factor l/y/2 and has an enhanced amplitude but that for 
the acoustic mode is almost dispensable. We plot the RSOI strength dependence of the frequency ratio for the optic 
mode in FigUJa) and for the acoustic mode in FigEJc). Both modes have the same dependence on a in the long 
wavelength limit q = as expressed by Eq. i|12|) . At finite wavevector, the frequency ratio shows strong deviation from 
that of q = for the optic mode but it is not the case for the acoustic mode. To see more clearly the dependence of 
lu/ujq on q, the DQW frequency ratio oj/ujq vs wavevector q for optic and acoustic modes are plotted in FigHJb) and 
(d) respectively with RSOI strength a = 1, 2.5, 3.5, 5, 7, 8, 9cvo for curves counted from top. The curves of the optic 
mode are similar to those of the SQW plasmon with minima near oj(q) = 2akp and maxima near 2akp. Because the 
DQW optic mode has a higher energy than the SQW mode, the minima shift to lower q values by a factor of \j\[2. 
The frequency ratio of the DQW acoustic plasmon mode, however, has a very weak dependence on wavevector. Inter- 
branch transitions can lower the optic plasmon frequency more than ten percent and the acoustic plasmon frequency 
less than three percent. 

The electron density dependence of the frequency ratio lu/loq is plotted for the optic mode in Fig. Ufa) and the 
acoustic mode in Fig. [Tfb) in a wavevector range up to q = 3 x 10 5 cm _1 . In the long- wavelength limit q = 0, it is 
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FIG. 7: lu/uo of (a) the optic and (b) the acoustic plasmons vs n e in a DQW with a/ao = 5. Results of q = (solid), 10 
cm -1 (dotted), 2 x 10 5 cm" 1 (dashed), and 3 x 10 5 cm -1 (dot-dashed) are shown. 

described by Eq.(ll) for both modes. Fine structure appears in finite q as inter-branch transitions are switched on. 
For the optic mode, with increasing wavevector, the ratio decreases at first and a dip appear about uj = 2akp which 
is n e = no for q = 10 5 cm _1 . Then the ratio increases from the low-density regime to the higher-density regime. At 
q = 3 it is higher than the value of q = over the whole regime. For the acoustic mode, on the other hand, the ratio 
lu/loq decreases and then increases with q in the low-density regime but keeps decreasing up to q = 3 in the higher 
density regime. The overall ratio at finite wavevectors is lower than that in the zero-wavevector limit. For n e > 3no, 
all the curves oj/ujq in Fig. J7J) are above 0.94 and converge to 1 with increasing n e . 

Note that our results are obtained in the RPA which may break down in systems of low electron density or systems 
at high temperature. The validity of the RPA can be justified by the parameter r s defined by n~ l — irr^a*^ with the 
effective Bohr radius a* B = 47re e oo ?i 2 /(m*e 2 ). For a 2DES of electron density no, which corresponds to r s ss 1.6 here, 
the RPA is still valid at low temperature. Experiments also indicate that the RPA describes the plasmon spectrum 
and the Coulomb drag of 2DEG's reasonably well at low temperature in GaAs/AlGaAs systems with electron density 
of order uq. [la. fl9L l21| For a system of lower electron density or at higher temperature, the local-field correction to 
the dielectric function and the plasmon spectrum of 2DEG's becomes more important. [l9l l20t l23l |24| 

IV. CONCLUDING REMARKS 

In the framework of the random phase approximation, we studied the dielectric function and the plasmon spectrum 
of single and coupled 2DEG's with RSOI in GalnAs-based semiconductor structures. In the SQW case, the RSOI 
lifts the spin degeneracy of electrons of finite wavevector but does not change the spatial wave functions when the 
spin orientation of a electron varies. Thus the bare Coulomb coupling between the electrons is not affected by the 
RSOI though the spin effect to this coupling should be considered explicitly in the presence of RSOI. The screened 
Coulomb coupling and the dielectric function, nevertheless, are modified by RSOI because the electron correlation in 
the system varies with the energy-band shape. As a result of the invariance of spatial wavefunction of electrons in 
different spin branches, the dielectric function describing the Coulomb coupling between electrons in different spin 
branches has the same form and the system is a single-component system. The deformation of the energy band from 
a parabola for each spin branch results in a negative term of intra-branch contributions to the dielectric function, 
which lowers the plasmon frequency. The frequency decrease is proportional to the square of the RSOI strength and 
the inverse of the electron density in the long wavelength limit. At finite wavevector of plasmon, the contribution of 
inter-branch transitions to the dielectric function is negative at energy lower than the spin split of the inner Fermi 
circle and positive at energy higher than the spin split of the outer Fermi circle. This leads to a shoulder in plasmon 
dispersion curve at the energy equal to the spin split of inner Fermi circle and a possible RSOI frequency enhancement 
for plasmon of higher energy dependent on the competition between the intra- and inter-branch contributions. 

In DQW systems, a similar influence of RSOI via intra-branch transitions is observed on both the optic and acoustic 
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modes. The effect of RSOI via inter-branch transitions is enhanced (diminished) for the optic (acoustic) mode. The 
inter-branch single-particle excitation spectrum covers part of the plasmon spectrum and leads to the Landau damping 
of plasmon. 
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